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Abstract — Most of the past research on economic 
replenishment quantity did not consider the trade credit period. 
In practice, however, companies in Taiwan have shown that 
taking the trade credit period into account when they are 
making inventory decisions decreases costs substantially. The 
aim of this paper was to develop an inventory model that 
incorporates a delay in payments for permissible shortages to 
minimize the total inventory-relevant cost. The optimal 
replenishment policy of the proposed model was identified by 
utilizing the global optimum conditions and non-constrained 
quadratic nonlinear programming. A numerical example was 
also generated. 

Index Terms — Optimal Replenishment Policy, Delay in 
Payment, Permissible Shortages. 

I. INTRODUCTION 

Inventory is one of the most important assessments when a 
firm attempts to improve the efficiency of the operations in its 
supply chain. Inventory control, which is linked to most 
aspects of production organization, can influence a company's 
business and has a critical role in the activities of operations 
management. Therefore, a company can take advantage of 
tracking stock by shortening the lead time for replenishment 
and estimating inventory cost as long as the company 
maintains a sturdy inventory management system. 

Goyal(1985) argued that IBM’s efforts in integrating their 
national spare part stock network, developing a new inventory 
system, and modifying their customer service quality have 
reduced their total inventory cost by more than $250,000,000 
and continuously saves the company up to $20,000,000 in 
annual expenses. These findings demonstrate the importance 
of inventory management in a company's business. 
High-working capital on inventory could result in 
expenditures such as holding costs and insurance, and the 
costs of other personnel matters might dramatically rise. Most 
manufacturers believe that flawlessly managing their 
inventory can not only reduce inventory cost but also increase 
their freedom when they use the capital in their companies. In 
contrast, shortages would occur if companies fail to control 
their inventory. Shortages would also lower their service 
levels, harm their goodwill and damage their competitiveness. 

Researchers like Wilson(1934) extended the original 

model and applied it to solve inventory problems in the real 
world. In addition, most Economic Production Quantity 
(EPQ) models incorporate the following assumptions: 

• Production rate must be greater than demand 

• There must be instantaneous replenishment 

In fact, the order quantity submitted by retailers usually 
limits the production problems that suppliers need to deal 


with. Therefore, it is necessary to relax the EPQ assumptions 
to create an inventory model that will find a better solution for 
the issues that suppliers face. Furthermore, demand and 
supply do not have to be balanced because shortages occur 
when an insufficient quantity of product is produced by 
suppliers. When the order quantity cannot be satisfied, the 
retailer either waits for replenishment or switches the order to 
other suppliers who offer the same products. 

Most of the traditional EPQ models assume that retailers 
pay for the goods immediately when the transaction is 
completed. In reality, few retailers pay for their goods 
immediately because most of them use trade credit. The use of 
trade credit violates the EPQ assumptions. Furthermore, to 
create a win-win situation, both the supplier and the retailer 
usually set up a reasonable credit transaction deadline by 
negotiating with each other. By setting a deadline, the supplier 
will accrue extra interest, and the retailer will not need to 
initiate loan financing. When delayed payment is provided by 
suppliers, the retailer also has the advantage of reducing 
interest yielded by the fund which uses in the payment of 
backlog. 

Based on the reasons we mentioned above, we constructed 
a theoretical model that incorporates a delay in payments 
under conditions that allowed shortages. We then 
demonstrated some important properties associated with 
implementing an optimal replenishment policy. After 
obtaining the model properties, we generated an example. In 
the last section of the paper, we provided our conclusions as 
well as suggestions for further research. 

II. Literature Review 

A. Credit Deadline for Payment Delays 

Goyal(1985) suggested an EOQ model with a fixed credit 
deadline for paying for goods as a promotion that could be 
provided by the supplier. If the retailer sells the product 
within the limit, he/she will be spared the interest generated 
by the funds that must be used to pay for the goods. That is a 
type of chance cost. If the pay date passes, he/she must pay for 
the products and the interest generated by the funds that have 
already been paid to the supplier for the rest of the products 
that are in stock. Many researchers have extended the above 
model in different ways. Chung(1998) investigated the 
interactions between credit transactions and the time value of 
currency, and he then employed a discounted cash flow 
(DCF) model to build a total variant cost function for 
inventory. In his research, an optimal replenishment policy 
was provided. Chung(2003) revised the model by Goyal(1985) 
and also found an optimal order quantity. Huang(2004, 2007) 
rewrote the assumptions made by Goyal(1985) and argued 
that the selling price and the purchase cost for the product 
might not be equal in real-world scenarios. Huang also 
considered the rate of replenishment to be finite and 
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developed his research in two phases: one where the supplier 
allows the retailer to pay for the products before the credit 
deadline, and another phase where the retailer also provides 
his/her customers with permission to delay payments for the 
sake of increasing the amount of product that is sold. 
Huang(2007) also discussed a possibility in which the 
supplier offers partial permission for delays in payment when 
the actual order quantity is less than a certain amount. He 
minimized the total cost to ensure that the inventory cycle and 
the order quantity were optimal. Liao(2007) investigated a 
non-instantaneous inventory model of deterioration for items 
that incorporated delays in payment. Tsao and Sheen(2007) 
found an optimal selling price and maximized the profit for 
deterioration items while including permission for delays in 
payment. Chang, Teng and Goyal(2008) reviewed the 
literature for inventory models published in the last twenty 
years and classified these models as non-deterioration items, 
deterioration items, models that allow shortages, order 
quantity and inflation. They also suggested further research 
directions for every model they reviewed. Ouyang, Teng, 
Goyal and Yang(2009) extended the work by Goyal(1985) by 
making a rule that the promotion cannot be offered if the order 
quantity submitted by the retailer is below a certain threshold. 

B. Permissible Shortages 

Shortages occur when production cannot satisfy the 
quantity of orders. If customers allow for delays in delivery 
within a reasonable deadline, it is a permissible shortage. 
According to Cohen(1977), complete backlogging means that 
customers will not switch to other suppliers or cancel the 
order because of loyalty to the product during shortages. In a 
competitive market, if the supplier cannot deliver products to 
the retailer and solve shortages, the retailer will likely switch 
his/her order to other suppliers, which results in a loss to the 
suppliers. Wee(1995, 1999) considered conditions that 
included partial backlogging and set up a constant between 0 
and 1 to be the backlogging rate. He also argued that waiting 
time was the most important factor that could affect whether 
the customers wait for replenishment or not. Philip(1974) 
argued that the deterioration rates of deterioration items 
should follow a three-parameter Weibull distribution. 
Papachristos and Skouri(2000, 2003) argued that the quantity 
that needs to be replenished is a decreasing exponential 
function of the time until product delivery. Jolai, 
Tavakkoli-Moghaddam, Rabbani and Sadoughian(2006) 
considered the factor of inflation and assumed that 
deterioration rates for deteriorated items should follow a 
two-parameter Weibull distribution to solve the EPQ model. 
Tsao and Sheen(2007) found an optimal selling price and a 
replenishment cycle with permissible delays in payment. They 
also sought to maximize the unit time profit. Dye(2007) 
determined a unitary selling price under an optimal 
replenishment policy. Dye(2007) also found that total profit 
was a concave function. Chung and Huang(2009) extended 
the work of Goyal(1985) by using dual variables and the year 
relevant cost to construct an inventory model. They also 
proved that the year relevant cost function was convex. A year 
later, Hu and Liu(2010) adopted the finite replenishment rate 
assumption and relaxed a condition in which purchase price 
and selling price must not be consistent as part of their effort 
to build upon the work of Chung and Huang(2009). 
Furthermore, they found an optimal inventory policy in the 
environment they had made. 


III. The Model 

A. Notation 

D Demand rate 

P Replenishment rate 

H Holding cost per unit of product in stock 

I max Maximum inventory level 

t The time point when the quantity of products 

in stock reaches I max 

t The time point when all of the products in 

stock are sold out 

t 3 The time point when the quantity of backlog 

reaches N 

T Inventory cycle 

lift) Total inventory in phase i and time point t, 
i= 1,2, 3, 4; t=t h t 2 ,t 3 ,T 
s Sale price 

Cl Shortage cost per unit of product 

N Maximum shortage level 

c Purchase price 

/ Inventory level at the time point of t — M 

I e Interest income for the unit fund 

4 Interest yielded by the unit fund, which is 

used to pay for the product in stock 
M Credit transaction deadline 

TVC(T) Total cost function 

B. Assumptions 

A3.1 Single product with non-deteriorating properties 

A3. 2 If both D and P are known and constant, it is reasonable 

to assume that P > D 

A3. 3 Purchase price and sale price of the product are both 

constant and s > c 

A3. 4 Instant replenishment 

A3. 5 I k ^I e 

A3. 6 When 0 — T—M, the retailer does not have to pay for the 
interest, 4, for the products in stock. Otherwise, when 

O^M^T, the retailer needs to pay the interest, I k . 

C. Model Formulation 

Based on the above notation and assumptions, the 
inventory model was developed and is presented in Figure 1 . 


Inventory level 



According to Figure 1 , we partitioned the model into the 
following phases. 

Phase 1: We start at t=0, which is the time point when the 
inventory level equals zero. As products are replenished at a 
rate P, the inventory level increases and reaches its maximum, 
Imax. Its increasing rate can be presented as follows: 
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dl { (t) 
dt 


= P-D, 0<t<t! 



Phase 2: Let I 2 ( t) be the inventory level when t^ [t } , t 2 /. The 
retailer stops purchasing any product, but market demand 
results in a gradual decrease in inventory level. At t= t 2 , the 
inventory level is reduced to zero. There is a variation in the 
inventory level during this period. 


dl 2 (t) 
dt 


= —D , L < t < t 2 



Then, the state of the inventory system enters Phase 3. 
Phase 3: Shortages happens in this period. The quantity of 
backlog reaches it maximum level, A, at t— t 3 Then, Phase 4 
begins. The rate of change of the inventory level during this 
period is: 


dl 2 (t) 
dt 


—D , t 2 < t < t3 



Phase 4: The retailer decides to replenish products to meet 
market demand. Because of the assumption P>D , the 
inventory level reaches zero when all backlogs are satisfied. 

dm =P _ D ;t3 < t < T ( 4 ) 

dt 


In the following content section, we investigate the optimal 
replenishment policy under the conditions 0 < M <t 2 and 
0<t 2 <M, respectively. 

The objective of the model is to minimize the total cost. 
Thus, we first defined the total cost as equal to the sum of the 
holding cost, shortage cost, the interest yielded by the funds 
used to pay for the product in stock and the negative interest 
income. By using (1) to (4) as well as some triangular 
properties, we found the following formulas for computing 
other cost functions and determining the properties of an 
optimal policy. 


t 


l 


pt 2 

p 




i 


max 



N=D { l ~j 

I„ =(h-M )D 


( 6 ) 

(7) 

(8) 
(9) 


According to Figure 1 and (5) to (9), the following cost 
functions can be derived. 

The holding cost during the period of 0 < t < t 2 is: 


h 


£' I { ( t)dt + 1 1 2 ( t)dt 


hDtl ( D \ 
1 

P 


V 


(10) 


) 


The shortage cost, which occurs during the period of 
t 2 < t < T, can be formulated as follows: 


a \j'rl I 3 (Odt-\- § I A (t}dt 


aD 

r i 

2 

l p ) 


](r - ) : 


(ii) 


Case 1 : 0 — M ^ t 2 

The interest yielded by the unit funds used to pay for the 

product in stock, I k , during the period M ^ t — 1 2 is shown in 
Figure 2. 


Inventory level 



Figure 2 The interest yielded by unit funds used to pay 
for the product in stock during the period 0<M<t 2 . 


Based on Figure 2, we constructed the interest yielded by the 
unit funds used to pay for the product in stock as follows: 

cL f 2 L ( t)dt = ( t - m) 2 (1 2) 

K JM z 2 

When 0 <M <t 2 , the interest income is constructed as 
written in formula (13). 

cl f 1 Pdt + si f 2 Ddt + / f (sD - cP)dt = I ((c + s)Du +-(sD- cP)( 1 -—)(T-t 2 ) 2 ) 

e Jo e Jo ’ e IP P (13) 

According to the aforementioned definition of total cost, 
we defined the total cost function as follows: 


hDt 2 

TVC l (t 2 ,T) = 2 


2 + 


^ (1 


cDI 


+ ^(t 2 - M y- 


I e ((c + s)Dt 2 + 


lf_ 

2 P 


(sD-cP)d 


D 7 

-)C T~t 2 f 


(14) 


Case 2: 0<t 2 <M 

Under this condition, retailers sell all of the products in 
stock before they have to pay for the products they purchased. 
Therefore, they do not have to pay interest. In other words, 
retailers are not affected by interest income during the periods 
(0, t 2 ) and (t 2 , M), as shown in Figure 3. 


Inventory level 

A 



Figure 3: The chance cost of interest income during 

0<t 2 <M 


For this condition, the total cost can be constructed as 
follows: 

hDtl 

TVC 2 {t 2 ,T) — 2 
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— (i - — )(r - 1 2 Y - 

2 P 


I e ((c + s)Dt 2 + ^ (?D- - c)( 1 - ~p)(T — t 2 ) 2 ) (15) 
D. Properties 

At the point when T = t 2 = M, if TVC\(M) = TVC 2 (M ), the 

cost function, TVC(T), is continuous. 

Based on the above assumption, we employed the optimal 
conditions for finding the replenishment time to minimize the 
total cost. 


dTVC l (t 2 , T) = ^ hDt ^ _ aD(T _ u ))(( _ 22) + c DI k (r 2 - M ) 

dt 2 - P 

- /, (Z>(c + s) - -T (Z ) 2 (sZ> - cP)(l -f )(T-, 2 „) (16) 


t 2 =(P 2 (aT + cl k M + I e (c + s + cDT)) - P(aDT + 7 e (DT(c + $)) + sl e DT) / 

(P 2 {a + h + cl t + cDI e ) - P((h + a)D + D 2 / e (c + s)) + jD 5 /, ) 

Let 

A = P 2 (a + h + cl k + cDI e ) - P((/z + a)D + D 2 I e (c + s)) + s£> 3 / e 

By substituting (17) with (14), taking a derivative with 
respect to the inventory cycle, T, and setting the resulting 
formula equal to zero, we can solve the equation and find the 

optimal inventory cycle, T x . 

a ' 7V f (r) = 4r (P\D 2 I e QiT(s + 2c) + csl k (T-M)-c 2 I k M-I e (c + s) 2 ) 
dT PA 

+ D(acl k (T - M) - al e (c + s) + lahT)) 

+ P(D 3 I ( si (c + s) - hT(c + 2s) - csL ( T-M )) - ah&T) 


e 

4 


+ shIDT) 


A = 


(cI k M + I e (c + s))P 


(18) 


(19) 


h(P — D) + cl k P 

To assure that 1\ is the global minimum, we took the 
derivative of (18) with respect to T again. 


d 2 TVC l (T) _ 1 


d z T PA 

\2 / t~\2 


( P 3 (a(h + cl, ) + cDI (h + cl, ) 


F (s - c)DI 

We then employed 1 k > again. It is clear that 

v 

A > 0 and (20) >0. □ 

By repeating the procedure for finding l\ we attempted 

to find the derivatives of TVC 2 ( t 2 , T ) and T 2 . 

By taking the first-order derivative of TVC 2 ( t 2 , T ) with 

respect to T, we obtain 

dTVC 2 (t 2 ,T ) 


D 


D 


dtr 


= hDt 2 (1 - — ) - aD(l )(T -t 2 ) 

- /, (Z>(c + 5) - ^ ( D 2 (sD - cPX 1 - “XT’ - f 2 » (21) 

When (21) equals zero, we obtain the optimal 
replenishment time after some rearranging. 


C = ( P 2 ((a + cDI )T + / (c + 5 )) - PT(aD+ ID 2 (c + s) + .v/DT) / 


((P - D)(P(a + /?)- / DOD - cP))) 


( 22 ) 


dTVC 2 (T) _ D_ 
dT _ 

v2t 2/ , x r rx2. 


Denote 5 = (P " D)(P(a + “ cP)) 

By taking the derivative of TVC 2 ( t 2 , T ) with respect to T 

and letting the resulting function be equal to zero, we can 
solve the equation to obtain the optimal inventory cycle. 

(— P~ ( cDI e 2 (s + c) + I e ( a(c + s) — chDT ) — ahT} 

+ P{sD l l l e {c + 5) - hDdlAc + s) - ahDT) + shD 3 TI e ) (23) 

ZZXlil ,24) 

/*(/> - D) 

Now, to examine the convexity of TVC 2 ( t 2 , T) , we took 

the second-order derivative of TVC 2 ( t 2 , T ) with respect to 
T. 

d TVC 2 (T ) ^hp_^_ p i^ a + cDI ^ + + DI (c + _ sTP I ) 

d~T PB 

Lemma 2: The conditions sufficient to ensure that 


r r~ i e 

7 2 — 


TVC 2 (t 2 , T 7 ) is convex 


are 


P> 




and 


/ > 


a(v - ch) 

cD(s + c) 

vD 


-P (D I (h(2c + s) + cl Ac + s)) + aD(2h + cl k ) Proof:LetP> — and P>0. 


2 , ( 20 ) 


+ P(L) I e ( h(c + 2s) + csl k ) + ahD " ) - shI e D" ) 

Lemma 1: The sufficient conditions under which 
TVC X ( t 2 ,T ) is a convex function are 

D ^ sD (s - c)DI 

P > and I, > - 


c s 

Proof: To show that TVC l (t 2 , T) is convex, we needed to 
prove that (20) is greater than zero. 

Because P > D , the condition of A > 0 is 

cP 2 I k -sD 2 I e (P-D)> 0 


To show that TVC 2 ( t 2 ,T ) is a convex function, we need 

to prove that (25) > 0, which requires a condition to be 
found such that 

- P 2 (a + cDI ) + hDP(a + DI (c + s)) - sD 3 I >0 Let 


P = 


sD 


and substitute it into the formula above. 


After rearranging, we have I > 


a(s - ch) 
cD{s + c) 


Let P = 


sD 


and substitute this into the formula 


Thus, if I > — — — — , the second-order derivative of 

cD{s + c) 


T (s - c)DI 

above. Then, we obtain I k > by utilizing (20). 


TVC 2 (t 2 , T) >0 


□ 


Denote 7J and T 2 as the extreme points of 
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TVC l (t 2 , T x ) and TVC 2 ( t 2 ,T 2 ), respectively. If Lemmas 
1 and 2 hold, we have the following solutions: 


TVC X (T) 


< 


< O , T 
= O , T 
> O , T e 


= (o,t; ) 

= T* 

(T * , oo) 


and 


tvc 2 (r) 


< 


<0, T 
= O , T 

> o ,r €= 


= (o,t; ) 

= T 
1 2 

(7" 2 * 9 00 ) 


(26) 

(27) 



IV. Numerical Example 


To obtain sufficient usage of their production capacity, 
suppliers will offer retailers a promotion for delaying 
payments until a credit transaction deadline, M. Therefore, the 
retailers can apply the results found in this paper to build their 
inventory model. Furthermore, an optimal replenishment 
policy that dramatically reduces total cost can also be 
obtained. Data were collected from some companies located 
in Taiwan and processed with the model. To mimic reality, 
the parameters used in the proposed model were assigned as 
follows: 


P = 50 

D-30 

H 

• 

o 

II 

II 

M =1.2 

a — 0.3 

s = 8 

o 

• 

o 

II 

04 

H 

• 

o 

II 

•-"H 



After performing the calculations, we obtained 


T* = 1.4625 . U= 2.25 w , M > u * a 

1 and 2 . When 2 , we found 

that the longer the optimal inventory cycle was, the greater the 
profit was that retailers could obtain. 


As a result, we plotted and TVC 2 versus T in Figures 


4 and 5, respectively. These two figures provide evidence that 
the aforementioned lemmas are valid. In addition, the delay in 
payments offered by certain suppliers would result in a 
negative inventory cost, which also leads to more products 
being purchased from the suppliers. 



Figure 4: The value of 7VC, versus T during 0 — M^t 2 . 


Figure 5: The value of TVC 2 versus T during 

0<r 2 <M, 


V. Conclusion 

The aim of this paper was to develop an inventory model 
that incorporates delays in payments for permissible 
shortages. The proposed model for an optimal replenishment 
policy was identified by utilizing the global optimal 
conditions for the Maple-based numerical analysis. We found 

that both TVC l (T) and TVC 2 (T) were convex functions. 

This result assures the existence of the polar point that is 
unique under certain conditions. Moreover, as shown in 
Figures 4 and 5, adopting the promotion offered by the 
suppliers would benefit the retailers’ businesses because the 
retailers would be more willing to purchase products from the 
suppliers. Therefore, a win-win scenario was achieved. 

Accordingly, in this research, we assumed that the product 
was not deteriorating and that no defective items were 
produced during the replenishment process. For future 
analyses, we will relax the assumptions by considering 
deteriorating items and assuming that the demand rate and 
replenishment rate follow certain probability distributions 
(e.g., a Weibull distribution or GAMMA distribution) to 
provide solutions for replenishment issues in the real world. 
In addition, taking the factors of inflection and time 
preference into account would also be worth investigating 
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